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For a particle traversing a bent crystal in the regime of volume reflection we evaluate the probabil- 
ity of interaction with atomic nuclei. Regardless of the continuous potential shape, this probability 
is found to differ from the corresponding value in an amorphous target by an amount proportional 
to the crystal bending radius, and the particle deflection angle. Based on this result, we evaluate the 
rate of inelastic nuclear interactions, and the flnal beam angular dispersion due to multiple Coulomb 
scattering. The theoretical predictions are compared with the experiments. The impact of multiple 
Coulomb scattering on the mean volume reflection angle is also discussed. 

PACS numbers: 61.85.+p, 29.27.-a, 45.10.-b 
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I. INTRODUCTION 

Deflection of fast eharged particle beams by bent crys- 
tals is a developing branch of accelerator technology, pur- 
suing the goal of steering ultra-high energy beams in a 
restricted laboratory space. Till the last decade, the con- 
ventionally exploited deflection mechanism in bent crys- 
tal tools was channeling but an alternative promis- 
ing technique based on volume reflection [2| has emerged 
recently, dealing with over-barrier particles, and offer- 
ing the benefit of large angular acceptance (equal to the 
crystal total bending angle) and almost 100% deflection 
efficiency. Some drawback of the volume refiection tech- 
nique is the smallness of the defiection angle, which is 
of the order of Lindhard's critical angle 9c = \/2Vq/E, 
with Vb the planar potential well depth, and E 3^ Vq 
the particle energy. This angle is independent of the 
crystal thickness, but a buildup of the deflection angle 
can nonetheless be acheived by transmitting the particle 
through a sequence of bent crystals [3|]; or arranging a 
composite volume reflection from several atomic planes 
in one crystal ^^l. An additional advantage of the vol- 
ume reflection mechanism is that it applies equally well 
to negative particles, for which channeling is spoiled by 
the unavoidable multiple Coulomb scattering on atomic 
nuclei residing at the bottom of the potential well. 

The origin of the volume reflection effect is essentially 
dynamical, due to non-linearity of the motion in the 
strong continuous potential of bent atomic planes. At 
that, incoherent scattering effects on individual atomic 
nuclei can be held small compared to the deflection angle 
mean value. That even admits using for volume refiection 
experiments the crystals of thickness 1-^2 mm, which by 
~ 10 times exceeds the volume reflection region extent, 
estimated as ^ ROc However, in other respects, 

for instance for evaluation of the deflected beam angular 
spread, the multiple Coulomb scattering can not be ig- 
nored. Worth mentioning also are the inelastic nuclear in- 
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teractions, which can be registered as multiple secondary 
particle signals in the beam loss monitors. Such events 
are relatively rare, and therefore present no danger for 
the primary beam propagation, but from their rate one 
can determine the robustness of the volume reflection in- 
side the crystal. In fact, since physically they occur due 
to interaction with the same nuclei as Coulomb scatter- 
ing, the rates of these processes must be closely related. 
Similar processes also include atomic K-shell ionization 
and the accompanying characteristic X-ray radiation oc- 
curring at fast charged particle passage close to an atomic 
nucleus. 

Outside the volume refiection area, we know that the 
particle motion becomes highly over-barrier, whereby the 
rate of fast particle scattering on atomic nuclei must 
approach that in an amorphous medium. Thus, in the 
thick-crystal limit, the number of nuclear interactions in 
the whole crystal roughly equals that in an amorphous 
target of same material and thickness. But that kind 
of approximation is unsuitable if the crystal thickness is 
comparable with the volume reflection area extent, or if 
one aims to exploit the nuclear interactions, for monitor- 
ing the in-depth particle dynamics. Since the number of 
nuclear interactions in an amorphous medium, is known 
precisely enough from various experiments, it may be 
neatly subtracted and in this way the difference be mea- 
sured. Yet the theory interpreting this difference needs 
to be developed. 

For a consistent treatment of incoherent multiple scat- 
tering at volume reflection, one needs to solve the kinetic 
equation in a non-uniform external fleld, but in general 
that does not seem feasible by analytic means. Still, there 
must exist a domain in which the Coulomb scattering is 
sufficiently small for perturbative account, and there are 
reasons to expect it to cover the practically interesting 
range of beam and crystal parameters where most of the 
presently available data belong. 

The aim of the present paper is to evaluate the rate 
of nuclear interactions adhering to the framework of the 
perturbative approach. Given the Id (radial) character 
of volume reflection dynamics in the pure continuous po- 
tential of a uniformly bent crystal, the corresponding par- 
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tide trajectories may be expressed analytically, and the 
probability of a nuclear interaction be calculated along 
the known path. In SecHIlwe implement this procedure, 
for particles of different charge sign, and different orien- 
tations of the crystal: (110) and (111)- In Sec. |lll] limi- 
tations of the perturbative approach are determined. In 
Sec. llVl we apply the results of Sec. |n]to the probability 
of inelastic nuclear interactions at specific experimental 
conditions, comparing the theory with the experiment. 
In Sec.|V]we extend the theory predictions to elastic scat- 
tering and address the issue of the volume-reflected beam 
angular divergence. The conclusions are given in Sec. lVIl 



II. INELASTIC NUCLEAR INTERACTION 
PROBABILITY UNDER NEGLECT OF 
COULOMB SCATTERING 

Volume reflection process assumes particle interaction 
with the bent crystal in a planar orientation. At that, 
the atomic density in each plane may be regarded as uni- 
form. Since all the nuclei are located in the planes, one 
particle crossing of an atomic plane may be viewed as an 
elementary act of nuclear interaction. If all the planes are 
packed them with the same density and equidistant [l3|, 
with the inter-planar distance d, the probability of any 
kind of nuclear interaction in one atomic plane crossed 
at a tangential angle 9 is 



Pi = "•atCTA- 



sin 9 



(1) 



where cfa is the corresponding cross-section on a single 
nucleus, and nat the atomic density in the crystal vol- 
ume. For elastic scattering one must employ the trans- 
port cross-section (a a = Ctr), while for inelastic interac- 
tions - the total inelastic cross-section on a Silicon nu- 
cleus {a A = Tinei)- Morc specifically, the inelastic in- 
teractions may be genuinely nuclear, involving multiple 
hadron production, or electromagnetic, e.g., ionization of 
the small-radius atomic K-shell|18| and the accompany- 
ing characteristic X-ray radiation. In the Glauber ap- 
proximation, all the high-energy nuclear and electromag- 
netic cross-sections are energy-independent, though one 
might take their weak energy dependence into account, 
if desired. 

In the case of a straight crystal, and for a highly over- 
barrier particle, when 9 is much greater than the critical 
value, 



> 



(2) 



and thus is subject to negligible variation within the crys- 
tal (var6' ^ 9c 9), summing up contributions ([ij for 
^'^"^ crossed planes would yield the total nuclear in- 
teraction probability: 



P = ni,t<JAL 



straight crystal, 
highly over-barrier motion 



(3) 



This value is independent of d and 9, and is equal to the 
corresponding probability in an amorphous medium - not 
surprisingly since the particle flow covers each nucleus 
with the same density, equal to that in the initial beam 
outside the crystal. In this sense, one can speak about 
an "amorphous orientation" of a perfect crystal even. By 
the same token, the latter notion applies in bent crystal 
regions where the particle motion is highly over-barrier. 

In the volume reflection case, however, condition ([2]) 
in the vicinity of a radial reflection point breaks down. 
In this region the plane crossing angle varies consider- 
ably along the particle path, hence it is to be evaluated 
accurately at each plane crossing. To this end, the par- 
ticle trajectory must be computed beyond the straight 
line approximation and, neglecting the multiple scatter- 
ing, we are going to approximate it by a trajectory in the 
pure continuous potential. If the crystal is bent uniformly 
(which is sufficiently credible nowadays), the trajectory 
description simplifies in cylindrical coordinates, with the 
bent planes corresponding to surfaces of constant radius 
relative to some axis far outside of the crystal. Thereat, 
the plane crossing angle sine entering Eq. ([T]) expresses 
simply as the time derivative of the particle radial coor- 
dinate: 



smt 



(4) 



Inserting Eq. (U) to Eq. ([1]) and summing over all the 
planes crossed by the particle, we obtain the total inelas- 
tic nuclear interaction probability in a bent crystal: 



P ~ natCTA 



(in a uniformly bent crystal) . 



(5) 

Since the particle motion straightens out away from 
the volume reflection area, there the nuclear interaction 
rate per unit length should approach that in an amor- 
phous medium. Hence, the difference between the num- 
ber of nuclear interactions in an oriented crystal and in 
an "unoriented" crystal may be expressed as 



AP = n^tfJAAL, 
where the isolated geometrical factor 

AL = lim 




L 



(6) 



(7) 



is expected to be finite and independent of L, represent- 
ing the excess (or deficit) of the target nuclear interaction 
range. 

In this section we will concentrate on evaluation of 
limit ([7]). To begin with, note that the trajectory in 
a centrally-symmetric continuous potential is symmetric 
with respect to the reflection point, enabling one to count 
the crossed planes beginning from the reflection point in 
only one direction and double the result: 



AL = 2 lim 



(8) 
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Here trcR is the distance from the reflection point to the 
crystal boundary - say, its front face, where the en- 
trance angle relative to the atomic planes is known 
(see Fig. [IJ. This quantity is expressible through the 
particle trajectory, too. In fact, at large rimax it is unam- 
biguously related with the total volume reflection angle 
^rofi, which can be evaluated in a model approach or 
measured. In the small-angle approximation, from Fig.[T] 
one infers ngl, 



9dcfl/2 = hm {00 - t,,ti/R). 



(9) 



Having traded trofi for Oq, the latter angle now is to be 
related with nmax- But this relation, in fact, appears to 
be trivial in the thick crystal limit whose condition is 



Vo 

m E 



("thick" -crystal limit). (10) 
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FIG. 1: Relation between the volume reflection angle, the 
angle of particle entrance to the crystal, the depth of the 
reflection point and the crystal bending radius. Dashed arcs 
- bent atomic planes. Sohd curve - schematic of the particle 
trajectory. 



Indeed, utilizing the transverse energy conservation law 
in a cylindrically symmetric continuous potential V{r)^ 



(11) 



we extract 



where 



(13) 



^ Ed E 
Rr ~ ~ cm — 

Wq 5 GeV 



is the Tsyganov's critical radius 1]. Solving Eq. p3|) for 
6*0, we get 



1 2d ^ 



O 



I R V 

2dnmax 2E 



(14) 



the inter-planar continuous potential. From Eq. (|lip we 
express 



Ei-V r 



E 



Rr 



(16) 



with E±_ = Ej^{0Q,b) depending on the particle initial 
conditions including its impact parameter b and the inci- 
dence angle with respect to the planes. Now, granted 
the periodocity of the intra-crystal continuous potential, 
values of V{r) are equal at atomic plane locations: 



V{rn) = V^lreat. plane = COUSt. 



(17) 



As for r in the centrifugal energy term in Eq. (|16p . its 
value at different atomic planes differs only by a multiple 
of d: 

rn = ''n,„a, + (f^max " n)d (n < rimax)- (18) 

With this in mind, we can write 



Here the interaction-dependent correction term asymp- 
totically vanishes as nmax oo, and may be omitted 
under the limit sign. 

Combining Eqs. ([51 [5]) and (IT^ . we cast AL in form 



/"max /(1/2R 

AL = ROdefi+^VWA lim V 



(15) 



A. Orientation (110) 



To proceed, we need to evaluate the terms of the se- 
quence of angles entering the denominator in Eq. (|15p . 
For a silicon crystal in planar orientation (110), r„ actu- 
ally appears to be a fairly simple function of the plane 
order number independently of the precise shape of 



where variable 



2-(r] + m), m = 0,l,2,.. 



(i^j_-V'(r„._))i? , r„ 



Ed 



(19) 



(20) 



accumulates all the dependence on the inital conditions. 
An important notice is that it belongs to an interval of 
unit length: 



^min V — ^min + 1, 



(21) 



in order to secure the relation min{r^} = max{f^_]^}. 

Substituting Tn^^^-m from Eq. (fT^ to Eq. (|T5l) . we 
cast it in form 



AL = i?0dcfl + V2RdC Q, r?) (Si (1 



10)), (22) 
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f(l/2,!?) 




FIG. 2: Blue curve - behavior of Hurwitz zeta function (|23|) . 
red curve - its approximation (|24p . black dotted curve - ap- 
proximation (|25p valid for 77 > 1. The integral from the func- 
tion over this interval equals zero (Eq. ((33|). 



or 

1 3 

f?min < ^71 < »?min + |, m ^ Vl + J- (27b) 

For positive particles, the probability of case (|27ap equals 
3/4, while that of case (j27bp is 1/4. For negative particles 
at i? > 4i?c the case (|27al) is realized with the unit prob- 
ability, because the higher potential barriers completely 
shadow the minor ones, even in spite of the centrifugal 
energy tilt. 

So far the formulation held for particles of any charge 
sign. But in what concerns the last unknown quantity 
Vmim the situation turns principally different for posi- 
tively and for negatively charged particles. Below we 
will analyze these two cases separately. 



where 



C. Positively charged particles 



= lim 



E 







(23) 



The latter function is categorized as Hurwitz (or gener- 
alized Riemann) zeta with the parameter equal ^ (for 



a general definition of C,{a,v) see 
purposes, it may be approximated by 



For all practical 



1 



1 



207 + 1 



7J + 1, V77, (24) 



obtained by application to the sum in (1231) (with the first, 
singular term being singled out) of the Euler-Maclaurin 
formula [loj . Furthermore, at 7/ > 1 it admits a simpler 
approximation 



1 



(77 > 1) 



(25) 



(see Fig.©. 



B. Orientation (111) 



In case of orientation (111), the continuous potential 
values at all the planes are still equal, but the intrepla- 
nar intervals assume alternating values d/A and 3(i/4. 
Thereat, the calculation principle remains the same, ex- 
cept that the summation over the atomic planes is carried 
out separately for odd and even numbers. The result then 
involves two different C functions: 



— ROdcR+y — 



(Si (111)), 
(26) 

where depending on which of the non-equivalent planes 
is encountered the last, 

?7min < ^yi < ?7min + |, »y2 = ?/l + | , (27a) 



In the case of positively charged particles (typically 
protons, which are most important for ultra-high energy 
accelerator applications), determination of 77min is par- 
ticularly simple. Consider, again, the case of orientation 
(110). At particle entrance to the reflection interval (see 
Fig. EK), the minimum of the kinetic energy is achieved 
when the particle passes the last potential barrier with a 
vanishing kinetic energy. But for positively charged par- 
ticles that barrier coincides with the atomic plane, the 
kinetic energy on which we want to know. So, at Tymin 
this energy merely tends to zero: 



Vn 



R 



2d' 



{rLJ 



0, 

(28) 

i.e., T] belongs to the interval < ?; < 1 (actually ex- 
hibited in Fig. [5]). The divergence of function ({1/2, rj) 
at the physical interval end-point 77 corresponds to 
a grazing crossing of the last atomic plane. The other 
end-point value equals C(l/2, 1) = C(l/2) ~ -1.46. Note 
that with the known lower limit (|28|) . 77 may even be 
found without the need to evaluate nmax, if one rewrites 
Eq. ^ as 



V = Ml 



iE± — V^|reat.plane)-R ''plane 

Ed d 



(29) 



with braces { }t signifying the operation of taking the 
fractional part, and rpiano being r at any atomic plane 
(without a difference under the { }f sign). 

Since at 77 < 1 typical values of function C are ^ 1, the 
ratio of the two terms in (12^ is of the order 



V2Rd 
R0C 



RcFir 



RVo 



jRc 
R 



(30) 



It is known that the quality volume refiection is possible 
provided R > AR^. (see Q), whereat ratio (PH)) appears 
to be < 1. But in actual practice that ratio may still be 
sizeable, especially at 77 — where C function blows up. 
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Veff(r) 




r, ,=V 2d/R , ;;„,„=1 



Veff(r) 



V|tc 



V|pla, 




'/max— ^min"'"! 



'■,..„„=-^2(l/|,„p-V|p|„e)/£ ^ !7m 



FIG. 3: Determination of the range of variation of parameter 
77 in crystal orientation (110); (a) - for positively charged 
particles; (b) - for negatively charged particles. For details 
see text. 



Then, it might seem that 77 must necessarily be specified 
issuing from the initial conditions. 

However, one should remember that in a real beam the 
initial conditions for the particle entering the crystal are 
not quite certain. The beam transverse dimensions are 
always much greater than inter-atomic distances in the 
crystal, therefore the particle impact parameters b have 
to be averaged over. Besides that, the indeterminancy 
SOq of the angles in the incident beam is usually large 
compared to d/L, and needs averaging in a vicinity of 
the mean value 9q. 

Examining Eq. ((29)) . we see that the expression un- 
der the fractional part sign contains contributions from 
the initial kinetic and potential energies in an additive 
way: E± = + V{b). Actually, fluctuations of both 
contributions are large: 



Ed 2 



R.9r 



S0n 



and 



R_ 

'Ed 



V 



RVo 



R 



Ed 4i?c 



> 1. 



(31) 



(32) 



Thus, at fluctuations of particle parameters in a real 
beam, 77 spans its unit definition interval many times, 
and each time its relation with 9q and b is approxi- 
mately linear, so in fact, 77 may be treated as a uniformly 
distributed random variable. Then, averaging over the 
beam is equivalent to unweighted averaging over r]. But 



such an average of function C entering (|22p gives zero due 
to identity 



(33) 



straightforwardly checkable from definition (1231) . That 
leads to a simple and model-independent relation: 



(AL) = R ( 



7dcfi/ 



(positively charged particles). 



(34) 

The same result (IM|) is obtained for orientation (111). 



D. Negatively charged particles 

The distinction of the negatively charged particle (typ- 
ically, TT~) case is that atomic plane positions do not co- 
incide with the tops of the potential barriers - thus, here 
Vmin > (see Fig.lSja). To determine the value of rymin, we 
need to know the particle kinetic energy at the last plane 
crossing, which in the present case equals the difference 
of potential energies between the atomic plane and the 
top of the preceding barrier: 



E . .2 
— mm r„ 
2 



— ^Iretop ^ ^IrGplane ~ Vq 



^ Rc 



(35) 

Here the last equality is exact for a parabolic inter-planar 
potential, while for a non-parabolic one it also turns to be 
exact in the limits R^ Rc and i? — > Rc, so heuristically 
we may expect it to be sufficiently accurate on the whole 
interval R > Rc- Therewith, we derive 



R 2 
77^,„ = _mm7^„_« — 



1 - 



Rc 



(36) 



In this expression, sa which for quality volume 
reflection is supposed to be > 1. From Fig. [2] we see that 
at 77 > 1 function C(l/2,77) is fairly smooth and may be 
linearized in 77 about the midpoint 77min + 1/2. Taylor- 
expanding Eq. (1^ . we get: 



1 



V - 7?min - 



(37) 



Averaging of the last term of Eq. (|57| over the unit 
interval of 77 gives zero. Substituting in the first term 
of (1271) 77,nin from Eq. and all that to Eq. dH]), we 

obtain the final expresson for the average nuclear range 
difference in the case of (110) orientation: 



(AL)=i?(0dcfl)-2i?0e(l-^ 



(38) 



(negatively charged particles, Si(llO)). 
At > 4i?c, when {9 

dcfi) ~ 9c, Eq. (j38p reduces to 
(AL) « -R9c. (39) 
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For orientation (111) in the negative particle case the 
calculation is more complicated. We will qoute the result 
under the condition R > 4i?c, retaining only the linear 
correction in Rc/R, which is relatively simple: 

(Ai) = R (^?defi) - 2R9c (^1 - (40) 
(negatively charged particles, Si(lll)). 



Here we let 0^ = ^JW^fE, and i?c ~ with Vl the 

larger well depth in a straight crystal (see [ij). 

Eqs. (j34[ I38p are our main results. Comparing them, 
we see that for positive particles the nuclear range excess 
is positive, while for negative particles it is negative (i. e. 
there is a deficit), being of the same order in magnitude. 
This situation is rather natural since positive particles 
are repelled from the atomic planes and cross them more 
tangentially, while negative particles are attracted, cross- 
ing the planes more steeply. From the practical side, it is 
worth noting that in order to detect the discussed effect, 
it might actually be easier to measure the difference not 
between the volume reflection case and the amorphous 
orientation, but between volume reflection cases for pos- 
itively and negatively charged particles in the same bent 
crystal - since thereat the signal is ^ 1.5 times larger. 



III. VIOLATION OF VOLUME REFLECTION 
DYNAMICS BY MULTIPLE SCATTERING 

The perturbative treatment of multiple scattering in 
the previous section permitted us to derive simple for- 
mulae, but their physical reliability yet depends on the 
crystal and the beam parameters. In this section, we will 
work out quantitative conditions at which the perturba- 
tive treatment of multiple scattering is justified, concen- 
trating on positive particle case (which is simpler and 
more important for applications). We will also touch 
upon the general trend in the f^dcfl dependence on R due 
to multiple scattering. 



A. Condition of coherent dynamics dominance 

The quantity representing the multiple scattering in- 
fluence on the volume reflection is the multiple scattering 
mean square angle (Tam within the volume reflection area 
extent ~ i?0dcfl- It competes with the dynamical an- 
gles, but it is a subtle issue to find among them the most 
sensitive quantity. For positively charged particles, most 
certainly, the smallest angle is that of the last atomic 
plane crossing, rn„,ax- Employing Eq. (|19p. the average 
'^"max i^^'Y be estimated as 



(41) 



Note that in contrast to Gc^ it does not depend on the 
particle energy, instead involving the crystal bending ra- 



dius. At particle energies and crystal radii suitable for 
volume refiection, {rura^.^} is smaller than 6c, as long as 



2Rc 
R 



< 1. 



(42) 



The angle of multiple scattering (in projection onto one 
relevant transverse direction perpendicular to the active 
family of atomic planes) is determined by the Highland- 
Lynch-Dahl formula [lj| [20j 



O-am(r) 



13.6 MeV 

E 



|- (^1 + 0.038 In |-), (43) 



with T the traversed material thickness, and Xq the 
material-dependent radiation length constant (for silicon 
Xq « 9.36 cm). The origin of the logarithm of T in 
Eq. is due to the Rutherford large-angle "tail" of 
multiple scattering, slightly affecting the Gaussianity of 
the profile. Formula (j43|) works with an accuracy of a 
few percent for high- and intermediate-Z substances. At 
T - 0.2 1 mm, i. e. 1 + 0.038 In -2^ ~ 0.8 ± 0.03, 
Eq. (|43|) may be used in the simplified form 



O'am(r) 



11 MeV 
E~ 



(44) 



(Note that the coefficient 11 MeV here is appreciably 
smaller than estimate ^f^^ = 2L2MoV _ 3 y^ y 

" Q \/2 



V2 



often used within the simplest leading logarithm approx- 
imation) . 

Inserting into Eq. (|44|) T ~ R6c, and dividing by 
Eq. (|4ip . we arrive at a requirement 



(TamiR6c) „ 11 MeV 



E 



Xod 



« 1, 



(45a) 



which may be viewed as a restriction on the crystal bend- 
ing radius: 

R <^ RinuitiE) (coherent dynamics dominance), 

(45b) 

if we introduce 



-Rmult 



E Xnd 



11 MeV V 9c 
(positively charged particles) 



E V" 



(We evaluate 0c = y/2Vo/E with Vq = 22.7 eV for 
Si(llO)). From Eq. ^ we note that at relativis- 

tic particle energies definitely i?muit >• 4i?c- Hence, there 
exists a range of crystal curvatures 4_Rc < R < i?muit cov- 
ering the practically most relevant cases R > ARc , and at 
the same time not subject to strong multiple scattering. 

It is also instructive to compare the restriction R < 
Rmuit with that stemming from the crystal finite thick- 
ness: 



R<^L/0c^ 100 m. 



(47) 
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m if 



Limitation (j45|) will be more restrictive than 

sis ^ ( i.75^Tev )^^^- Hence, at not very high energies 
and not very thin crystals, such as they are nowadays, 
boundary effects on the volume reflection are inessential. 
In principle, the boundary condition effect is to reduce 
(^dcfl), but the multiple scattering can manage that alone 
- as we will discuss below. 



(''defl)/Slim 

1.0 ' 



0.8 
0.6 
0.4 
0.2 




400 GeV 



- I 
13 GeV 



B. Suppression of the mean angle of volume 
reflection by multiple scattering 

The effect of multiple scattering on (^dcfi), which is 
needed for input in Eqs. (|341 138p . in the first order must 
vanish, due to the symmetry of the incoherent scatter- 
ing in the scattering angle sign. Hence, (^dofl) is also a 
convenient variable for representing the onset of multiple 
scattering non-linearity. 

In general, the strong multiple scattering manifests it- 
self as follows. Once condition (j45b|) is violated, the 
multiple scattering can affect the particle trajectory in 
the last volume reflection interval. At still larger R, it 
may affect the trajectory even before the last interval 
is reached. But since the non-zero mean defiection angle 
effect receives significant contribution from the last inter- 
val, if the particle does not reach it in due way, the vol- 
ume reflection may be effectively terminated. This must 
be deflnitely true for those particles which scatter out- 
wards in r, while those scattered inwards may still have 
a chance to volume reflect later. Thus, at i? > i?muit 
the value of (^dcfi) is from general reasons expected to 
decrease. Eventually, the particle net deflection ought to 
vanish as the crystal straightens out. But the rate of the 
decrease can not be assessed without a proper calcula- 
tion, or reference to experimental data. 

On the other hand, it is well known that in the region 
R < ^muit function (^defl) {R) grows with the increase of 
R. There exists a deflnite theoretical prediction for the 
mean angle of volume reflection in a pure continuous po- 
tential (applying harmonic approximation to the latter, 
which works well in (110) orientation for silicon) f7l|. For 
positive particles the formula reads [2T| 



7dcfl/harm 



71im 



1 - 



d 



2' 



(48) 



FIG. 4: Dependence of the mean volume reflection angle for 
positively charghed particles on the crystal curvature. Blue 
points - experimental data for protons at 400 GeV llj. Solid 
curve - theoretical prediction (|48|l taking into account only 
dynamics in the pure continuous potential. Red points - mea- 
surement for positively and negatively charged particles at 13 
GeV [l^]. Dashed line - interpolation. 



For positively charged particles and orientation (110), in 
the range R > i?muit there is only one experimental point 
- at i5 = 13 GeV and R — 2.4 m, which corresponds to 
R ~ 3.6i?„iuit- On the other hand, at i? < i?nmit a rel- 
atively detailed experiment was performed at E = 400 
GeV |11| . To link these two experiments, one may sug- 
gest that at i? > i?muit ^ Rc, scale Rc drops out, and so 
(^'dofl) /6'iiin becomes a function of the ratio R/Rmuit{E) 
only. Then, in the plot of (6'dcfl) /6'iim vs. R/Rmuit data 
at 13 GeV may be used to continue the dependence at 
400 GeV (see Fig. [4]). Its behavior, slowly decreasing 
with R beyond i?muit, is close to that obtained in sim- 
ulation @ at a single energy 400 GeV [2^. Thus, our 
expression (|46|) for i?muit(-E') for positive particles is not 
unreasonable, and so is the scaling law hypothesis. 

For negatively charged particles, it is unobvious 
whether formula (|46)) is applicable as well. Even if not, 
there must exist a similar expression for i?muit for nega- 
tive particles, monotonously increasing with E. However, 
there are experimental indications (l2l . [l3| that the vol- 
ume reflection pattern at the same energy for positively 
and negatively charged particles strongly differs. For neg- 
ative particles with the decrease of energy, and hence 
increase of R/Rmuit, the mean deflection angle rapidly 
diminishes. We will not speculate about this behavior 
until more detailed evidence arrives. 



(4i?c < R < i?muit, positive particles) . (49) 

Since function (^dcfi) (R) decreases to the both sides 
of i?muit, it must achieve a maximum somewhere around 
^^muit - a prominent feature to look for in computer sim- 
ulation and experiments. 

C. Scaling law at large R 

At the present state of affairs, the data on volume re- 
flection are not very abundant, particularly at i? > i?muit- 



IV. COMPARISON WITH EXPERIMENT FOR 
INELASTIC NUCLEAR INTERACTION 
PROBABILITY 

We are flnally in a position to test predictions of 
Sec. [H] against the available experimental data. The 
most straightforward check is supposed to be with the 
results of experiments on inelastic nuclear scattering. At 
present, the only such experiment is that with 400 GeV 
protons using a L = 2 mm thick silicon crystal at a single 
value of the crystal bending radius i? = 10 m Q. 

The results of this experiment are displayed in Fig. [S] 
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interval SOq ~ 2 /irad. From Eq. (|3ip we get 



67^ = ^00(5^0 - ^SOo ~ 10. 
a 2a 



(52) 



Apparently, this number is still much greater than the 
unit interval of variable 77 definition, hence with the vari- 
ation of 9o in that range, 77 actually scans its defini- 
tion interval several times, i. e. cutting angles down 
to ~ 1 /irad still can have no significant impact on AL. 
So, we can only attribute the small-cutting-amgle fluctu- 
ation of (AL) to enhanced experimental errors reflecting 
the difficulty of achieving such small cutting angles. 



FIG. 5: (adapted from 6]). Dependencies of the inelastic 
nuclear interaction probability of 400 GeV protons in the R = 
10 m crystal on the cutting angle of the incident beam: (1) 
- for 'amorphous orientation', (2) - for the case of volume 
reflection. 

(adapted from j[6||). In fact, the inelastic nuclear in- 
teraction rate was measured for varying cutting angles 
SOq. When the cutting angle was sufficiently large (which 
ought to correspond to perfect averaging over b or E±), 
the measured relative difference was about constant, on 
the level 

« (5 ± 2)%. (50) 

Our prediction, using the experimentally determined 
mean value (0defl)exp = 13.35 ^rad at the given curvature 
l/R = O.lm"^ amounts tofH] 

m^RV'.f...^.^^ (51) 

The theoretical accuracy of prediction (|5ip can be esti- 
mated as 0[R/ Ri-anXt) (see Eq. ([45]) ). Substituting here 
i? = 10 m, and 

^muit(400GeV) « 50 m, 

we infer J?/i?muit ~ 1/5. That is commensurable with the 
relative difference between our theory and the experiment 
(Eqs. (15(11) and ([CTt ). Needless to say, the approximate 
agreement of one number does not conclusively validate 
the theory. Measurements at different crystal thicknesses 
and bending radii would be worthwhile. 

A curious feature in Fig. [5] is that at small cutting 
angles (< 2 /irad) the difference between the rate of in- 
elastic nuclear interactions at volume reflection and in an 
amorphous case seems to depart from a constant and ac- 
tually vanish, although experimental errors in this region 
are too high for an unambiguous conclusion. In principle, 
some sensitivity to 59^ could emerge due to imperfect av- 
eraging and to the impact of the second term in Eq. (1^^ . 
To check this possibility, let us estimate the range of vari- 
ation of the argument rj under the variation of 6q in an 



V. ANGULAR SPREAD DUE TO ELASTIC 
MULTIPLE SCATTERING 

Registration of inelastic nuclear scattering discussed in 
the previous Section requires dedicated instrumentation 
like beam loss monitors. But even in a minimal beam 
deflection setup, the nuclear interactions shall manifest 
themselves through an angular broadening of the final 
beam. A complication here arises because the broaden- 
ing receives an additional contribution from the impact 
parameter dependence of the defiection angle, even in a 
pure continuous potential. In fact, the latter contribu- 
tion is anisotropic, but the beam spread transverse to 
the direction of deflection is rarely measured, so in the 
published experimental data on the beam dispersion in 
the direction of deflection they contribute together. As 
usual, the main source of broadening is away from the 
volume reflection region, but we suppose it to be sub- 
tractable. The problem, again, is to explain the differ- 
ence from the amorphous orientation, and the difference 
between positive and negative particles. 

The proper subtraction is possible if the kinetics of the 
particle passage through the crystal is decomposed into 
three distinct stages: pure incoherent multiple scattering 
upstream the volume reflection region (where the beam 
acquires Gaussian shape), pure dynamical broadening in 
the volume reflection region, and pure incoherent mul- 
tiple scattering downstream of it. In terms of the cor- 
responding angular distribution functions that expresses 
as 

_(9-92)^ _ ef 

de J ' V2^a2 J d{e,-e,) ^a,' 

(53) 

where we assumed the upstream and downstream inco- 
herent scattering to be purely Gaussian, and dwcoh is the 
angular distribution function in a pure continuous poten- 
tial (slightly averaged over the particle incidence angles). 
The small portion of multiple Coulomb scattering within 
the volume reflection region may be included either in 
the upstream or downstream scattering piece, as long as 
it is small, and thus additive. Moreover, if in (|53l) we 
change the integration variables, one integration can be 
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taken, with the result 
dwcohia) 



dw 



—T- = da 



da 



(54) 



depending only on the sum cr^ + o-| . Obviously, the latter 
sum must be equated to (t^^{L + AL), and thereby the 
precise positions of the boundaries separating the differ- 
ent kinetic regions prove to be inessential. 

Examining Eq. ([5l)) . one should realize that for real 
crystals, usually, the width of dwcoh/da is smaller than 
crfjjj, whereby the resulting angular distribution becomes 
close to a Gaussian, anyway. So, it is described essentially 
in terms of two moments; 



and 



de{e - {9)y 



dw 



(55) 



(56) 



The mean value {9) only receives contribution from 
dWcoh/dO: 



(0) 



doe- 



,dWr 



dO 



(57) 



(as was implied in Sec. while when we evaluate 
from Eq. (|54p , the coherent and incoherent contributions 
to it appear to be just additive: 



AR (^dofi)) 



'^coh' 



with 



coh 



d9 io-{d)y 



dw, 



coh 



de 



(58) 



(59) 



Angular distribution dwcoiJdO was evaluated in Q 
in the model of harmonic continuous potential between 
(110) silicon crystallographic planes. It has some dif- 
ferences, for positive and negative particles, as does the 
nuclear interaction rate calculated in Sec. |lTl Let us be- 
gin with the case of positive particles, to which most of 
the data refer. 



A. Positively charged particles 

For positive particles, at i? > ARc the coherent part of 
the angular distribution looks as (see Eq. (72) of [7]) [24] 



dwc 



dO 



Rdc Q ( TTd 

nd " \2Rec 



\o-{o) 



(60) 



Q being the Heavyside unit step function, and {9} being 
given by Eq. (|48|) . The corresponding (Tcoh ensues as 



fcoh 



(61) 



notably being inversely proportional to the crystal bend- 
ing radius. For a realistic continuous potential the nu- 
merical coefficient in (|6ip may slightly differ, but that is 
not crucial for the following estimates. 

Measurements of total for 400 GeV protons inter- 
acting with a (110) silicon crystal were carried out in 
experiment ,11]. There, in order to get access to the 
intrinsic volume reflection angular divergence Ucoh, the 
difference 



(62) 



was evaluated. From Eq. ((58|) we see, however, that it 
differs from pure ficoh: 



'coh 



7^ CTcoh- 



,(i?(0dcfl)), (63) 



Assuming condition (j45a|) to hold, we may insert 
explicit theoretical expressions (|44l) . ((6T|) . (|48|) into 
Eq. (|63l) . which leads to a non-scaling i?-dependence of 
the measured quantity fiv.r.: 



T2elR^^~\ E ) 

{R < -Rmult)- 



Xn 



(64) 



The most characteristic feature of function ((64|) is the 
existence of a minimum. The minimum location is found 
by equating to zero the derivative of the radicand with 
respect to R: 



R-AE) 



1 



-Xod^ 



E 



\ 2/3 

llMeV/ 



38 



E y/' 

-:z m. 

GeVy 

(65) 

The physical meaning of i?, is not vastly different from 
that of i?inuit ~ it marks the scale of R where the multiple 
scattering compares with coherent deflection angles, with 
the proviso that i?muit is derived from generic reasoning 
in terms of the particle trajectory, while i?* - in terms 
of specific contributions to the beam broadening in the 
crystal. However, the actual expressions for i?muit and 
Rf: differ; moreover, their ratio 



mult 



TTllMevj \d) \l2~\5OMev) 

(66) 

even depends on the particle energy, albeit pretty weakly. 
At ultra-relativistic energies ratio is > 1.5, justifying 
the use of Eq. (|64|) . but up to LHC energies it does not 
exceed 3. The minimal value of a^.i. 



Cmin = CTv.r.iR = R*) - 

sets the scale of a^.r. at i? > i?* 
function (j64[) varies slowly. 



ikcvy/^ 



E 



(67) 



because beyond i?* 



10 



The available data at E = 400 GeV (presented in 
Fig. [S]) do not reach beyond i?muit, but do reach well 
beyond i?,, which according to Eq. (|65l) equals 



i?*(400 GeV) w 16 m. 



(68) 



Around the data show a flattening of the R- 
dependence, but the point R — 35.71 m closest to i?muit 
seems to resume the decrease of CTv.r., which is rather un- 
expected. It must be noted that at this value of R the 
final beam shape qualitatively changes - there develops 
a thrust in the elsewhere Gaussian profile, extending to 
the side of the crystal bending, though not as far as to 
the channeling angle L/2R (see in ll| Fig. 3, left panel). 
That component is not associated with volume reflec- 
tion, and is conventionally attributed to volume capture, 
but still remains a poorly understood phenomenon (see, 
e.g., [1, [l^). Besides particle capture to the channeling 
regime, one can in principle imagine other mechanisms 
of particle drag to the side of the crystal bending, e. g., 
being due to detention of a fraction of protons on the top 
of a curved potential barrier, containing an area of strong 
multiple scattering. In either case, it is not very surpris- 
ing that an additional component arises as R approaches 

■Rmult ■ 

Concerning the observed reduction of a^.r. compared 
to the theoretical prediction ([M]) at i? = 35.71 m, one 
might say that in view of the final beam profile non- 
Gaussianity, as well as the ambiguity of the separation 
of volume reflected and volume captured fractions, our 
description neglecting the volume capture is invalidated 
as a whole, at least for description of CTv.!-.- On the other 
hand, the volume captured fraction comprises only 6% 
of the particles, and the volume reflected fraction is still 
Gaussian, whereas the discrepancy with the theory by a 
factor of ^ 2 appears to be exceedingly large. 

In the search of an explanation to the encountered dis- 
crepancy, we must pay attention to the fact that at eval- 
uation of a difference between two close quantities the 
experimental errors enhance. Evaluating propagation of 
errors for a^.r. defined by Eq. ([63]) . we strictly derive 



a/(ct ± (5(t)2 - (cTam ± <5f7am)^ = CTv.r. ± (5(Tv.r. , (69) 



with 



(70a) 



-V(W+W^ (if 



(70b) 

Now, for the point at i? = 35.71 m the prefactor 
Cam/o'v.r. in Eq. (j70b|) amounts ^ 5, but the experimen- 
tal error for (Tv.r. is quoted to be as small as Sa^.i-. — 
VO.lP + 0.092/irad — 0.14/irad, being comparable to 
(^Cam — 0.11/irad. Thus, probably, the error propaga- 
tion factor of 5 yet needs to be added for this point, and 
for consistency - for the preceding two points as well[25|. 




R. 20 



R(m) 



FIG. 6: Subtracted final beam angular width vs. the crystal 
bending radius, for E = 400 GeV protons in a L = 2 mm 
silicon crystal. Solid curve - theoretical prediction (Eq. (I64[l ). 
Dashed curve - with experimental (6^deH)- Dotted curve - pure 
(Tcoh- Blue points - experimental data [ill]. Gray error bars - 
evaluated by Eq. (fTO)) . 



The error bars modified in this way are displayed in Fig.[S] 
by gray color. The theoretical curve, in turn, may shift 
somewhat lower (see the dashed curve in Fig. [5]) if instead 
of ^lini one uses the empirically determined (^dofl) , as in 
Fig. m Therewith, the discrepancy between experiment 
and the theory appears to be just within one standard 
deviation. Formally, it may even comply with cTcoh as 
well (the dotted line in Fig. (|6])), but in the light of the 
foregoing analysis, that seems unlikely. 



B. Negatively charged particles 

Experimental data for negative particles at present are 
too scarce to allow phenomenological investigations of 
tj^ J. , so we will describe for them the situation in general 
terms. 

First of all, for negatively charged particles the ex- 
pression for crcoh(-R) somewhat differs, though its R- 
dependence remains close to 1/R, up to logarithmic fac- 
tors. Secondly, in this case we have A(7a„i cx (AL) < 0. 
Therefore, the expression for ^. for negative particles 
is similar to the radicand of Eq. (|64p , only with a nega- 
tive coefficient at the second term. That implies that for 
negative particles ^ changes sign and becomes nega- 
tive for sufficiently large R. That is the salient feature of 
the final beam angular distribution for negative particles, 
which would be interesting to verify experimentally. 

The second remark is that since (Tcoh for positively and 
for negatively charged particles differ, in general it is not 
as straightforward to compare the broadenings for posi- 
tive and negative particles, as that was the case with the 
rate of inelastic nuclear interactions. However, in the re- 
gion R> Rc where (Tcoh gets relatively small, that must 
already be possible. The simplest way of pinning down 
(Tcoh, though, is to measure both angular beam broaden- 
ing components perpendicular and parallel to the family 
of the active atomic planes. 
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VI. SUMMARY AND CONCLUSIONS 

In the present article we were concerned with the prob- 
lem of incoherent nuclear scattering at volume reflection. 
This problem belongs to the category of combined poten- 
tial and stochastic motion, which is complicated even in 
the radial Id case, hampering advances without a numer- 
ical simulation. However, it turned out that small (per- 
turbatively tractable) multiple scattering limit covers the 
practically most important region Rc R < RmuitiE), 
with Rmuit{E) defined by Eq. (|46l) . There, by analytic 
means we attained a simple relation that the difference 
between the probability of nuclear interactions of a pro- 
ton in a bent crystal and in an amorphous target is pro- 
portional to the product of the beam mean deflection an- 
gle and the crystal bending radius (Eqs. (O, ([M)) . ([5^ 1. 
That relation exploits only the local character of nuclear 
interactions and the periodicity of the atomic planes, and 
does not resort to any parameterization for the continu- 
ous potential, such as the parabolic approximation used 
in our earlier treatment of the volume reflection 0. 

Our perturbative prediction for positive particles was 
confronted with the experimental data @, [ll|, both for 
the rate of inelastic nuclear interactions (Fig. [S|) and for 
the volume-reflected beam broadening (Fig. [5]) . The the- 



oretical predictions are in an encouraging agreement with 
the data, though some unclear phenomena still remain on 
edges of measurement intervals. 

The utility of the predicted correction is twofold. 
Firstly, it may be just a noticeable correction which needs 
account. With the trend to decrease L in order to reduce 
the background from useless multiple Coulomb scatter- 
ing, the relative significance of the crystal curvature de- 
pendent nuclear scattering must increase. Secondly, upon 
accumulation of our experience on behavior of volume re- 
fiection at different experimental conditions it may grow 
into a method of monitoring the state of the volume re- 
fiection dynamics inside the crystal. The guiding princi- 
ple here must be, basically, that the greater the nuclear 
scattering difference, the more robust the volume reflec- 
tion dynamics. 

In conclusion, let us point out that besides the domain 
R < i?muit which matters most for the accelerator ap- 
plications, it is also of scientific interest to understand 
the particle passage process under conditions R > i?muit, 
where volume refiection merges with the volume capture. 
In Sees, mil IVl we tried to give some qualitative insight 
into these issues. The results obtained in the present 
article may also serve as a guide for more fundamental 
computer simulations. 
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